A Discrete Nonlinear Model with Substrate Feedback 
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We consider a prototypical model in which a nonlinear field (continuum or discrete) evolves on 
a flexible substrate which feeds back to the evolution of the main field. We identify the underly- 
ing physics and potential applications of such a model and examine its simplest one-dimensional 
Hamiltonian form, which turns out to be a modified Frenkel-Kontorova model coupled to an extra 
linear equation. We find static kink solutions and study their stability, and then examine moving 
kinks (the continuum limit of the model is studied too). We observe how the substrate effectively 
renormalizes properties of the kinks. In particular, a nontrivial finding is that branches of stable 
and unstable kink solutions may be extended beyond a critical point at which an effective intersite 
coupling vanishes; passing this critical point does not destabilize the kink. Kink-antikink collisions 
are also studied, demonstrating alternation between merger and transmission cases. 
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I. INTRODUCTION 

Soft condensed matter systems, such as vesicles, mi- 
crotubules and membranes, have recently attracted a lot 
of attention in both biological and industrial applications 
[1-4]. More generally, physical systems at the nanoscale 
including nanotubes and electronic and photonic waveg- 
uide structures [5,6] have nontrivial geometry and are 
influenced by substrate effects. These wide classes of 
problems, many of which are inherently nonlinear, raise 
the question of the interplay between nonlinearity and an 
adaptive substrate, including, in particular, a possibility 
of developing curvature in the substrate due to feedback 
forces from the nonlinear system. Some of these systems 
such as, for instance, the DNA double strand [7,8], are 
also inherently discrete. 

In this situation, the derivation of evolutionary equa- 
tions for intrinsic fields in a nonlinear system interfacing 
with a flexible substrate should take into regard the local 
dynamics and the feedback of the substrate. A context 
where these phenomena can easily manifest themselves 
is elasticity. Usually, models of particles connected by 
springs, such as the Frenkel-Kontorova (FK) model [9], 
assume that the springs are strictly linear. However, in 
reality the spring constant may depend on the spring's 
stretch. In this case, we are meaning not a straightfor- 
ward generalization of the FK model that includes the 
spring's anharmonicity, but a different system, in which 
the spring constant obeys its own dynamical equation, 
see below. 

There is an increasing body of literature dealing with 
the interplay of nonlinearity, discreteness, and a sub- 
strate subject to curvature. Usually, even if the substrate 
is curved, its geometry is assumed to be fixed, see, e.g., 
Refs. [10]. However, for many applications, ranging from 
condensed matter to optics to biophysics, it is quite rele- 



vant to examine a prototypical model that admits a pos- 
sibility of a variable substrate, which is affected by the 
primary fleld(s) and feeds back into its (their) dynamics. 

In this work, we introduce a simple one-dimensional 
model which captures some of the physically essential fea- 
tures of the interplay of a nonlinear system with a flexi- 
ble substrate. A natural setting for the model is a lattice 
with on-site nonlinearity, or its continuum counterpart 
(see below), hence the first ingredient of the system is 
the FK model (alias the discrete sine-Gordon equation 
[11]): 



CA2U„ — sinu,i 



(1) 



where C is the coupling constant (the spring constant, if 
we consider adjacent sites as being coupled by springs), 
A2W„ = {un+i + Un-i — 2un) is the discrete Laplacian, 
while the last term in Eq. (1) is the on-site nonlinearity. 
The overdot denotes the temporal derivative, while n in- 
dexes the lattice sites. The FK model has a plethora of 
physical realizations, the simplest one being the Scott's 
model, i.e., a chain of pendula suspended on a torque- 
elastic thread [12]. 

The springs can be made nonlinear by assuming that 
the constant C in Eq. (1) is replaced by a site-dependent 



one, C„ = Co - 



where Co is the constant mean value. 



while Vn is the variation of the spring constant due to the 
variation of the displacements Un and m„_i at the sites 
that the spring couples. The simplest nontrivial possi- 
bility, provided that we aim to produce a Hamiltonian 
system, is to assume that Vn responds to a change in 
displacements according to the following equation: 



k {Un 



l)] 



(2) 



where a is an intrinsic stiffness of the spring (see below), 
and k accounts for its susceptibility to the stretch. We 



1 



use the eonvcntion (sec also a sketch of the system in Fig. 
1) that the n-th spring eonneets the sites n — 1 and n. To 
maintain the Hamiltonian character and self-consistency 
of the model, the equation for the field u„ should be 
modified [cf. Eq. (1)] to read: 



= CqA 



2Ur. 



+ k{Vn 



Vn+1) 



(3) 



In fact, the extra terms added to Eq. (3) can be easily un- 
derstood in terms of the mechanical model: a difference 
in the value of the string constant produces a difference in 
the elastic forces even if there is no change in the length 
of the springs. Equations (2) and (3) are derived from 
the Hamiltonian 



1 (-2 



I 



-vl + (1 - cos(u„)) 



l)^ - kVn {Un - Un-l) 



We note that a somewhat similar (but still significantly 
different) two-component model was proposed to describe 
the coupling between protons and heavy ions in Ref. [13] 
(see also references therein). A single-component model 
with nonlinear springs was examined (mainly in terms of 
its thermodynamic properties) in Ref. [14]. 



The continuum limit (CL) of Eqs. (4) and (5) is 

Utt=Uxx-smu- kva;+r{vux)x, (6) 
vtt = -av + kux + ^ul, (7) 



where the constants k and r were appropriately rescaled, 
and the subscripts stand for partial derivatives. Setting 
r = in Eqs. (6) and (7) yields the CL of Eqs. (2) 
and (3). Note that, even in the latter case, the Lorentz 
invariance of the CL equations is broken (in comparison 
with the sine-Gordon equation, which is the CL of the 
usual FK model). An exact solution for the static kink 
(the topological soliton in the M-subsystem) can be easily 
found (for the case r = 0; a is then rescaled to be 1), 



M = 4 tan ^ 
v = 2 



exp 



X — Xo 



k ^ / X — Xo 
sech ' 



Vl - fc2 VVl - k"^ 



(8) 
(9) 



but exact solutions for moving kinks are not available 
even in the CL limit. 

Also worth noting is the CL for the linear spectrum of 
extended waves, which can be obtained by substituting 
u ~ Aexp(i{Kx — fit)) and v ~ B exp{i{Kx — ^It)) in 
the linearized version of Eqs. (6) and (7). For r = and 
with the normalization a = 1, we obtain 
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FIG. 1. Sketch of the model. Vn is the deviation of the n-th 
spring constant from Co due to the discrete gradient of the 
displacement Un- 

While this work is focused on the simplest model based 
on Eqs. (2) and (3), which includes only the linear cou- 
pling between the fields m„ and we note that a rele- 
vant generalization may involve quadratic couplings: 

iln = CoA2U„ - sin(M„) + k{vn - v„+i) 

+rVn{Un-l - Un) + rVn+l{Un+l - U„) , (4) 



aVn - k{Un- Un-l) 



- [Un - Un-l) 



(5) 



The origin of the quadratic terms is evident: they stem 
from second-order effects combining the change in the 
string length and elasticity. 



which consists of two phonon bands, one with ^J\ -\- \k\ < 
f2 < oo, and the other one with a/1 - jfcj <n < 1. No- 
tice the presence of a finite gap between these bands, 
provided that k ^ 0. In the case jfcj = I, the lower 
branch of the linearized spectrum is acoustic, while the 
upper one is always of the optical (alias, plasmon) type. 

Returning to the discrete case, there is no explicit so- 
lution for the static case; however, after the substitution 
of Vn = {k/a) {un — Un-l), which follows from Eq. (3), 
the steady-state equation for u„ becomes 



(Co - kya)A2U„ 



smur, 



(10) 



and hence k effectively renormalizes the lattice coupling. 
In fact that is what can be observed from Eqs. (8) and 
(9), as well as from Fig. 2: the effect of k is to decrease 
the width of the kink (by a factor approximately equal 
to jfcj, for k small). As is well known [11,15], a narrower 
kink is less mobile. A numerical solution of the static 
discrete equation (10) yields, similar to the case of the 
regular FK model [11,15], two types of steady-state dis- 
crete kinks, namely a stable inter-site centered one (an 
example for Co = 2fc = 1 is shown in Fig. 2) and an 
unstable site-centered one. 
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FIG. 2. A stable static intersite-centered kink («„ and Vn 
components are shown by circles in the top left and top right 
plots, respectively) for Co = 1 and k = 0.5; solid lines con- 
necting circles are guides to the eye. For comparison, the 
fields arc also shown for the case fc = 0, which corresponds 
to the usual FK model (stars connected by solid lines). The 
bottom subplot shows the spectral plane (Ar, A,) of linear sta- 
bility eigenvalues A found from the linearization of Eqs. (2) 
and (3) around the kink (for k = 0.5; the subscripts i and 
r refer to the real and imaginary parts of the eigenvalues). 
The absence of eigenvalues with nonzero real part implies the 
stability of the kink. Here, a = 1. 

We have examined the variation of the statie kink so- 
lutions and of their linear stability eigenvalues (for small 
perturbations) as a function of |fc|. Typical examples 
are shown in Fig. 3 for the stable and unstable kinks for 
Co = a = 1. A surprising finding in both cases is that the 
solutions could be continued for some range beyond the 
point |fc| = 1, which is effectively equivalent to the zero- 
coupling (anti-contrnwitm) limit, see Eq. (10). The shape 
of the solutions found for |fc| > 1 is non-monotonic (but 
they remain stable). These solution branches terminate 
at \k\ PS 1.04 for the stable solution, and at « 1.07 for 
the unstable site-centered one. 
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FIG. 3. The top subplot shows the intersite-centered solu- 
tion and its stability (top left and right panels) for \k\ = 1.04 
(just prior to the branch termination). The bottom panel 
shows the two stable internal-mode eigenvalues of the inter- 
site-centered kink (circles connected by the solid line) vs. 
k. The bottom subplots show the same for the unstable 
site-centered solution (the two top panels are for k = 1.07, 
i.e., also just prior to the termination of the corresponding 
branch, and the stars in the bottom subplot show the un- 
stable eigenvalue). Note that the profiles of both kinks are 
non-monotonic. 

The above-mentioned expectation that the kink be- 
comes less mobile as |fc| increases is verified dynamically 
by taking an initial kink boosted to the speed c = 0.2. 
Solving Eqs. (6) and (7) (with r = 0) with this initial 
condition for different values of A;, we notice that, as is 
shown in Fig. 4, the kink travels farther for k = than 
it does for k ^ 0, despite the fact that the initial velocity 
is the same in all the cases. In other words, the effective 
velocity which the kink demonstrates differs from that 
which was lent to it initially, and it also depends on |fc|. 




FIG. 4. Trajectories of a continuum-limit kink, initially 
boosted to the velocity c = 0.2, obtained from the numeri- 
cal integration of Eqs. (2) and (3) at different values of fc. 
fe = 0, fc = —0.05, fc = —0.15, and fc = —0.25 correspond, 
respectively, to the solid, dashed, dash-dotted, and dotted 
lines. The trajectories are shown only for late stage of the 
time evolution. 

Finally, we have also examined kink-antikink collisions 
in the system of Eqs. (2) and (3). An example of 
this type is shown in Fig. 5. In this case, Cq = 16/9 
{h = 0.75), and the initial speeds of the kink and an- 
tikink are ±0.1. The figure indicates the possibility 
of existence of parameter windows for which the kinks 
eventually escape each other's attraction (upon multiple 
bounces). Such windows lie between intervals of merger 
(kink-antikink annihilation into a breather). This feature 
is reminiscent of well-known findings of Ref. [16] for kink- 
antikink collisions in continuum (nonintegrable) models. 
However, unlike what was found in that work, here the 
windows are in terms of the parameter k, observed at a 
fixed value of the collision velocity. 

An example of such a window (a chacteristic case of 
which is shown in Fig. 5) occurs for 0.176 < k < 0.179, 
even though in the FK model (e.g., for fc = 0, a case 
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also shown in Fig. 5), for v = 0.1, only annihilation 
is possible. Notice that for this value (0.1) of the ini- 
tial speed (and within our resolution steps of 0.001 in 
k), no additional windows were identified. However, by 
considering different values of the initial speed (such as 
e.g., V = 0.125), we have verified that typically multiple 
such windows can occur. The dynamical behavior of the 
present model, which is essentially richer than in its FK 
counterpart, is provided for by the fact that the kink's 
internal mode is affected by variations of the coupling to 
the substrate (see the discussion above). 



alization of the model to two dimensions, which may be 
very relevant for physical applications. Furthermore, dis- 
sipativc versions of the model can be of direct relevance 
to chemical and biophysical applications. Detailed con- 
sideration of higher-order feedback effects, collisions be- 
tween kinks, and long-lived breather-like states in this 
model would also be of interest. 

This research was supported by the US Department of 
Energy under the contract W-7405-ENG-36. 




FIG. 5. Kink-antikink collisions for Co = 16/9 {h = 0.75) 
and initial speeds ±0.1 are shown in terms of trajecto- 
ries of centers of the colliding kinks. The solid, dashed 
and dash-dotted linos correspond, respectively, to fe = 0.0, 
k = 0.178 and k = 0.2. 

The case of a = in Eq. (2) has its own inter- 
est, as in this case w„ may be regarded as an acous- 
tic phonon branch (the one without a gap in its dis- 
persion relation), which interacts with the optical (plas- 
mon) branch represented by the field u„. Equations (2) 
and (3) with a = furnish an elementary example of 
this type of interaction. In this case, the equations have 
solely uniform zero- velocity solutions. Nontrivial solu- 
tions may only exist with a finite speed c. In particu- 
lar, the model's CL gives rise, in this case, to a double 
(47r) kink (not shown here), which is produced by the 
equation (1 — c^)w^5 — sin(u) -|- (fc^/c^)w^ = 0, where 
^ = X — ci. However, this simplest version of the model 
with a = is fiawed, as its Hamiltonian is unbounded 
from below, hence the model is not a well-posed one. It 
can be amended to fix this problem, keeping the optical- 
phonon and acoustic-phonon character of the fields Un 
and Vn, but detailed consideration of this issue is beyond 
the scope of the present work. 

To conclude, in this work we have proposed a simple, 
prototypical model, that may be developed to describe 
dynamics of kink-shaped excitations in more complex 
systems. A noteworthy extension would be the gener- 
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